Critical point in the strong field magnetoresistance 
of a normal conductor/perfect insulator/perfect conductor composite 
with a random columnar microstructure 
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A recently developed self-consistent effective medium approximation, for composites with a colum- 
nar microstructure, is applied to such a three-constituent mixture of isotropic normal conductor, 
perfect insulator, and perfect conductor, where a strong magnetic field B is present in the plane 
perpendicular to the columnar axis. When the insulating and perfectly conducting constituents 
do not percolate in that plane, the microstructure-induced in-plane magnetoresistance is found to 
saturate for large B, if the volume fraction of the perfect conductor ps is greater than that of the 
perfect insulator pi. By contrast, if ps < pi, that magnetoresistance keeps increasing as B 2 without 
ever saturating. This abrupt change in the macroscopic response, which occurs when ps = pi, is 
a critical point, with the associated critical exponents and scaling behavior that are characteris- 
tic of such points. The physical reasons for the singular behavior of the macroscopic response are 
discussed. A new type of percolation process is apparently involved in this phenomenon. 
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Theoretical and experimental studies, performed since 
1993, have shown that the macroscopic electrical re- 
sponse of three-dimensional (3D) composites with a two- 
dimensional (2D) or columnar microstructure can exhibit 
surprising forms of behavior when subject to a strong 
magnetic field B. For example, when a periodic array of 
parallel cylindrical holes is etched into an otherwise ho- 
mogeneous free-electron-conductor host, the system ex- 
hibits a strong dependence of the macroscopic or bulk 
effective resistivity on both the magnitude and the direc- 
tion of B, when B and the average current density (J) 
both lie in the plane perpendicular to the columnar axis. 
This occurs whenever the Hall-to-Ohmic resistivity ratio 
of the conducting host, H = PHail/pohmic = mI b I = w cT, 
(here fj, is the Hall mobility, lo c is the cyclotron frequency, 
and t is the conductivity relaxation time) is greater than 
1, and appears even if the array has a high point symme- 
try, e.g., square or triangular .Erl3 A strongly anisotropic 
magnetoresistance was also found in calculations on a pe- 
riodic columnar array of perfectly conducting inclusions 
embedded in a similar host.B More recently, composites 
with a disordered columnar microstructure were consid- 
ered, using an appropriate modification of the Brugge- 
man self-consistent-effective-medium-approximation.a In 
such a system, if the pure constituents have an isotropic 
transport behavior and the microstructure is isotropic in 
the plane perpendicular to the columnar axis, and if B 
lies in that plane, then both the longitudinal and the 
in-plane-transverse components of the bulk effective re- 



sistivity tensor p e 
(i.e., (J) 



denoted by p|, e) (i.e., (J) || (E) || B) 
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E) and they both lie in that plane 
but are perpendicular to B) respectively, are independent 
of the direction of B or (J) in that plane. In that study, 
it was found that for a two-constituent metal/insulator 



[M 1 1) mixture of this kind, both of these resistivity com- 
ponents increase as H 2 for \H\ ^>> 1, without ever reach- 
ing saturation. This behavior is not confined to the case 
where M represents a free-electron-conductor. All that 
is required is that its transport behavior be isotropic, 
and that the Hall-to-transverse-Ohmic resistivity ratio 
H have a magnitude that is much greater than 1. An 
experimental consequence of this nonsaturating behavior 
would be that the bulk effective magnetoresistivities of 
such an M/I mixture would continue to increase as B 2 
when B is large enough so that the Ohmic resistivities 
of the conducting constituent are saturated but its Hall 
resistivity continues to increase as B. By contrast, in a 
normal conductor/perfect conductor (M/S) mixture of 
this kind, both of those resistivity components saturate 
at finite values when \H\ ^> lJa 

Here we consider the macroscopic response of a 
three-constituent columnar composite of isotropic nor- 
mal conductor, perfect insulator, and perfect conductor 
(M/I/S), where the 2D microstructure in the plane per- 
pendicular to the columnar axis is again random, and a 
strong magnetic field B is applied in that plane. 

We use the "columnar unambiguous self consistent ef- 
fective medium approximation (CUSEMA)" , which was 
developed in Ref. m In this approximation, the self con- 
sistency requirement is that the in-plane components of 
the extra electric field E produced by isolated circular- 
cylindrical inclusions of the different constituents, em- 
bedded in the fictitious uniform effective medium, vanish 
on average. A similar requirement is not imposed upon 
the columnar component of E, which would be unmea- 
surable in a thin film realization of such microstructures. 
Also, the columnar component of (J) is always assumed 
to vanish. These requirements lead to the following self 
consistencey equations (we take x to be the columnar 
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axis) 



(I - Pe/ Pine) ■ 7inc(/0inc, Pe) ' Pt 



(1) 



where / is the unit tensor, ( ) indicates an average over 
the different types of inclusions, while the subscript {yz} 
indicates that only the y, z components of the 3x3 ten- 
sor appearing in the square brackets [ ] are included in 
this calculation, i.e., the 2x2 matrix in the lower right 
corner of [(I - p e /p inc ) ■ 7i nc (p inC! p e ) ■ p e \. The 3x3 ten- 
sor 7inc(Pinc! Pe) gives the uniform local electric field E int , 
which appears inside an isolated circular-cylindrical in- 
clusion, with resistivity tensor pi nc , inside the otherwise 
uniform effective medium host p e , when a uniform elec- 
tric field Eo is applied at large distances 



Ei, 



lb 



En 



(2) 



This tensor must be calculated from pi nc and p e — this was 
done for the relevant types of inclusions in the Appendix 
ofRef.|. 

It is important to note that the self consistency re- 
quirements used to derive CUSEMA differ from the re- 
quirements which are used to set up the conventional 
two-dimensional, Bruggeman-type, self-consistent effec- 
tive medium approximation (SEMA) in the y, z plane of 
this system (see Ref. [| for a discussion of SEMA in the 
presence of a magnetic field, when the conductivity ten- 
sors are nonscalar tensors) . Consequently, it should come 
as no surprise if the results of those approximations are 
sometimes quite different, as shown in Ref. ||. The great 
advantage of using CUSEMA is that the exact relations 
which p e must satisfy as a result of the duality symmetry 
are satisfied automatically, whereas they_can be seriously 
violated in the two-dimensional SEMA.Q 

For the resistivity tensors of the effective medium host 
and the three different types of inclusions we write 
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pi = pil, pi > pm; ps = psl, ps<Pa/- 



(3) 
(4) 



These forms mean that a magnetic field B is applied 
along the z axis, which is perpendicular to the colum- 
nar axis x. 

The form assumed for pm means that the M con- 
stituent exhibits isotropic transport behavior even when 
B is very strong, and does not rule out a dependence of its 
transverse and longitudinal Ohmic resistivities, pm and 
vpM, and also of its Hall coefficient Hp M /\B\, upon the 
magnitude of B. However, it does rule out the possibil- 
ity of clean single crystal inclusions of a transition metal 
like Copper, where the resistivities become very sensi- 
tive to the crystal orientation with respect to B and J at 
strong fields and low temperatures due to the existence of 



open orbits on its Fermi surface, although polycrystalline 
inclusions of such a metal might still be allowed. By con- 
trast, the form assumed for p e means that we expect the 
bulk effective transport behavior to be at least mildly 
anisotropic as a result of the columnar microstructure, 
i.e., in general we expect to have a e ^= j e . 

The forms assumed for pi and ps include the case 
where pi — oo and ps = 0. In fact, it is only for those 
extreme values that a mathematical singularity will actu- 
ally be found to appear in the bulk effective macroscopic 
response, justifying the term "critical point" . 

The same resistivity scale pu is used in both pm and 
p ei whose in-plane components will obviously have simi- 
lar magnitudes when neither the / constituent nor the S 
constituent percolates in that plane, i.e., when their vol- 
ume fractions pj, ps satisfy pj < 1/2, ps < 1/2. How- 
ever, because there are perfectly conducting inclusions 
that span the system from end to end along the colum- 
nar axis x, therefore the components of p e which involve 
that axis will all vanish. More precisely, we will have 
a e = 0(ps/pm), Pe = 0(ps/pm) as long as the volume 
fraction p$ of the S constituent is nonzero. 

The 7inc(Pinc, Pe) matrices for the [three types of inclu- 
sions which must be considered area (we omit the sub- 
script e from a e , (3 e , j e , and A e in these expressions, in 
order to save space) 



linc(PM,Pe) = 
( 1 





(-) (-) 








1/2 



V 



(•5) 



where we omitted the term f3' 2 /a, which is 0{ps / Pm) <C 
1, when it appeared alongside 7, which is 0(1), 



/ 1 



7inc(/3/,Pe) 



m l/2 










(6) 



i+m 1/2 y 

where we also discarded 0(pm / Pi) <C 1 terms, and 
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i+m l/2 



7mc(As,Pe) = 



/ 1 

' -I 1+(^) 1/2 



V 



PS 





PS 





(7) 



where we kept some 0(ps / Pm) -C 1 terms, because those 
terms are subsequently multiplied by an 0(pm/ps) 3> 1 
term — see Eq. (0). 

There will be two coupled nontrivial self-consistency 
equations, resulting from the yy and zz components of 
Eq. (Q) (the yz and zy components of that equation are 
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satisfied automatically due to the 2D isotropy of the mi- 
crostructure). Thus, we get two coupled algebraic equa- 
tions for the two unknown quantities A e = p^ / Pm and 

_ / 

7e = Pi VPm: 



\ 1/2 

0= I [(l-2p s )(l + H 2 )-(l-p I ) le +p I \ e ] 

- Ps (l + H 2 ) + (2 PI - l) 7e + (1 - Ps )(l + H 2 ), 

e 

(8) 



1/2 



VPS ~ lePl 



Kl ~PS) - A e(l - Pi)' 



(9) 



These equations can be transformed into a single poly- 
nomial equation for, say, A e . That would be an 8th order 
equation, which is quite complicated and which I have 
not been able to factorize algebraically. On the other 
hand, if we are interested only in the asymptotic behav- 
ior of A e and 7e when \H\ 3> 1, then an explicit solution 
can be obtained using asymptotic analysis. The results 
are 



' H^T^'-TT **Pl>Ps, \(pi-Ps)H\»1 7 

'Ptr^r t°*pi<ps, |(pj-p s )ff|»i, 



V^l^l^f forp/<Ps, \(pi~Ps)H\ « 1, 



7e 
Xe 



^f) for 7e »l, 



(i^f) 2 for 7e = 0(l). 



(10) 
(11) 



These results can be recast with the help of two scaling 
functions 



7e 



PI -ps 



-F(Z), A e 



Pi - ps 



-G(Z), (12) 



where the scaling variable is Z = \H\(pi — ps)/V^, and 
the scaling functions have the following limiting forms 



for Z > 1, 



F{Z) 



( 1-pj z 2 

pi l-2p r 

-^f (1 -2p s ) forZ«-l, 



(13) 



G(Z) 



i-p lz 
pi 

Pi z 2 

l-pi l-?Pl 



for \Z\ < 1, 
for Z > 1, 



< 



1— ps 



(l-2p s ) forZ«-l, 



(14) 



1— Pi 



z 



for \Z\ < 1. 



Both of these functions have the qualitative form shown 
in Fig. 0. 

Clearly, \H\ = oo, ps = pi defines a line of critical 
points of the macroscopic magnetotransport of such sys- 
tems: For ps > pi, both p^ and ~p~^ saturate when 
\H | — > oo, whereas for ps < pi they both keep increasing 
as H 2 for jiJj ^> 1. As ps — > pi from below, the coef- 
ficients of the H 2 terms tend to zero as pi — ps, while 
if ps — > from above the saturated values of both pft 

and p^f 1 diverge as l/(ps — Pi)- When ps = Pi, then 
these resistivity components keep increasing as \H\ for 
\H\ 3> 1. In the vicinity of the critical line, where both 
l/\H\ and \ps — pi\ are very small, it is easy to see that 
F(Z) S G(Z) and 7e = A e . 

Because these results were obtained within the frame- 
work of CUSEMA, some of these behaviors are not ex- 
pected to be correct in detail. We do expect that even 
a more accurate calculation of asymptotic behavior (i.e., 
^> 1) will exhibit saturated behavior for pj < ps, 
nonsaturating cx H 2 behavior for pj > ps, and nonsatu- 
rating oc \H \ behavior for pi = ps- However, we expect 
that more accurate calculations will show that the critical 
behavior as ps — > pi is not characterized by the simple 
forms l/(ps — Pi) or pi — ps which were obtained here, 
but rather by some noninteger values of the critical ex- 
ponents. Such calculations are now in progress. 

The behavior found in these calculations can be un- 
derstood qualitatively by recalling that the in-plane y, z 
components of &m are 

Pm\ ■ 



(15) 



For \H\ ^> 1, this represents a very anisotropic 2D con- 
ductor in the y, z plane, with pm&Mzz = 1/f = 0(1) 
and pM&Myy — l/H 2 C 1. In order to get from end 
to end of the sample, the in-plane electric current must 
make its way between different S inclusions by flowing 
through the M host. If there are many more S inclu- 
sions than / inclusions, then straight line S-to-S trajec- 
tories can easily be found that are parallel to z, allowing 
the current to flow through the M constituent only in 
the z direction. Therefore the macroscopic conductivity 
will depend only on <tmzz-> and the macroscopic resis- 
tivity will saturate when \H\ 3> 1. In the opposite case, 
when there are many more / inclusions than S inclusions, 
the current will often not be able to flow even between 
neighboring S inclusions only along z, but will have to 
have a nonzero y component in the M constituent. This 
is the low conductivity direction, therefore the macro- 
scopic conductivity will now be determined primarily by 
ouyy — 1/ (pmH 2 ). Threfore the macroscopic resistivity 
will not saturate, but keep increasing as H 2 forever. 

As the relative proportion of the I and S constituents is 
varied, there will be a transition from saturating to non- 
saturating cx H 2 behavior, which will have to be abrupt, 
i.e., it will be a singular or critical point of the macro- 
scopic response. This will be reflected by a similarly 
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abrupt change in the detailed local current distribution, 
which will have only z-parallel flow lines in the M con- 
stituent for values of pi below the transition point value, 
but will also have nonzero values of J y in that constituent 
when pj is above that value. Such a transition constitutes 
a new type of percolation phenomenon, which we believe 
deserves further study, both theoretical and experimen- 
tal. 

Experimental study of the critical point we have dis- 
covered could be done using a doped semiconductor film 
as the M constituent, with a random collection of etched 
perpendicular holes as the / constituent, and a random 
collection of perpendicular columnar inclusions, made of 
a high conductivity normal metal, playing the role of the 
S constituent. We would like to note that extremely low 
temperatures or very clean single crystals would not be 
required in order to observe this critical point. What 
would be necessary is a large contrast at each stage of 
the following chain of inequalities 

ps < PM < H 2 p M < pi. (16) 

If Si-doped GaAs is used as the M host, with a negative 
charge carrier density of 1.6 x 10 18 cm -3 and a mobility 
(i = 2500 cm 2 /Vs at a temperature of 90 K, as in the 
experiment described in Ref. |2j then a magnetic field of 
40 Tesla would result in H = —10. Such a material would 
have an Ohmic resistivity of 1.6 x 10 -3 f2cm, about 1000 
times greater than that of Copper. Thus, using Copper 
for the S inclusions and etched holes for the / inclusions, 
there should be no difficulty in satisfying all the above 
inequalities. 




FIG. 1. Qualitative shape of the two scaling functions F(Z) 
and G(Z)— see Eqs. @ and @. 
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